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The asymptotic behavior of the vorticity for the steady incompressible Navier- 
Stokes equations in a two-dimensional exterior domain is described in the case where 
the velocity at infinity is nonzero. It is well known that the asymptotic behavior 
of the velocity field is given by the fundamental solution of the Oseen system which 
is the linearization of the Navier-Stokes equation around Uqo- The vorticity has the 
property of decaying algebraically inside a parabolic region called the wake and ex¬ 
ponentially outside. The previously proven asymptotic expansions of the vorticity 
are relevant only inside the wake because everywhere else the remainder is larger 
than the asymptotic term. Here we present an asymptotic expansion that removes 
this weakness. Surprisingly, the found asymptotic term is not given by the Oseen lin¬ 
earization and has a power of decay that depends on the data. This strange behavior 
is specific to the two dimensional problem and is not present in three dimensions. 
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1 Introduction 


The stationary flow of an incompressible fluid past a body B is described by the Navier-Stokes 
equations, 

Ati — Vp = u • Vii , V • u = 0 , 

u\qq = u* , lim w = Uoo , 

|a?|^oo 

in the domain H \ H, where f is the source force, u^o 7 ^ 0 the velocity at infinity and u* 
is any boundary condition with no net flux. 



(lb) 


We assume that the body B is an open bounded domain with smooth boundary. In view of the 
symmetries of the equation, we assume without lost of generality that u^o = 2ei and 0 E H. 

This system has been subject to many investigations, see Galdi (2011, Chapter XII) for a 
complete statement of the main results known for this problem. Leray (1933) has shown the 
existence of weak solutions, but with the procedure he used, he was unable to verify that u 


1 


tends to at large distances. Gilbarg &: Weinberger (1974, 1978) have shown that any Leray 
solution u either converges at large distances to some constant vector uq, or the average of u of 
over circles in the L^-norm diverges as the size of the circle grows. Later on, Amick (1988) proved 
that if / = 0 and u* = 0, then u G L°° and therefore u converges to a constant uq at infinity. 
However, the question if uq = u^o is still open in general. In case Moo 7^ 0; Finn &: Smith (1967); 
Galdi (1993, 2004) used the Oseen approximation and a fixed point technique to prove existence 
and uniqueness of solutions to (1) for small data. The asymptotic structure of the solutions was 
presented by Babenko (1970) who shows in particular that velocity behaves at infinity like the 
Oseen fundamental solution. The asymptotic expansion of the velocity was also given under 
more general assumptions by Galdi & Sohr (1995); Sazonov (1999). The asymptotic behavior 
of the vorticity was first given by Babenko (1970, Theorem 8.1) only in the wake, and then by 
Glark (1971, Theorem 3.5’). These two results are relevant only in the wake region, i.e. for 
— < 1, because otherwise, the remainder is larger than the asymptotic term which is given 

by the Oseen linearization. In fact, we prove that the true asymptote which is also valid outside 
the wake region is not given by the Oseen linearization. Under smallness conditions, we show 
that the asymptote of the vorticity is given in polar coordinates (r, 6) by 

= 7-"^(l-cos0)+Bsin0 

for all e G (0,1), where A, B G M depend linearly on the net force F and /r is a 27r-periodic 
function depending on / and u*. Surprisingly the power of decay of the asymptote depends on 
the net force F and in particular this contradicts the statement of Theorem XII.8.4 in Galdi 
(2011) for any solution with F' / 0. 

Notation For x,xq G we use the following notation 

r = \x\, ro = |a;o|, ri = \x-XQ\, 

9 = Za;, 6q = ZiCo , 6i = A{x — xq) , 

where Zx denotes the angle 6 G (—7r,7r] such that x = \x\ (cos 0,sin 0). For a positive function 
re : H ^ M, we write A{x,xq) = 0{w{x)), if for xq in a bounded domain, there exists C > 0 
such that for all a; G H, 

|A(a;,a;o)| < C'|rt;(a;)| . 


rl/2 


o 


„l/2+e 


r(l—COS 0 ) 


2 Asymptote for the linear problem 


It is well-known that the problem (1) is related to the Oseen system which is the linearization 
of (1) around u = Woo = 2ei, 


Au — Vp — 2diu = / , V • u = 0 , 

U\g^ = U* — Woo ; li™ W = 0 , 

|a;|—>-oo 

The fundamental solution of the Oseen system is given by 


with 


/ di'tp - G d2'4! \ 

02 ^' -<91 " 0 j ’ 


il) 


H + G 
2 


= — log r , 


e = -VF, 



27r 


Ko{r) . 


( 2 ) 
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Denoting by T{u,p) the stress tensor, 

T(ii,p) = Vu + (Vu)^ — pi , 
the Green identity for the Oseen operator is 


I (V • T{u,p) — 2diu) • V — I (V • T(t;, q) + 2div) • u 
Jo, J S7 

= / {v • T{u,p) — u • T(t), q) —2u • V ei) • n . 
Jan 

Therefore, the solution of the Oseen system is given by 

/ E(a; — •)f — / E(a; — •) {T{u,p) — 2u® ei) + u • T(E, w){x — •) 

Jn Jan 


uix] = 


with 


/ u • T(E, w){x — •) • n = / [n • VE(a; — •) • u + u • VE(a; — •) • n — w{x — •)u • n] , 
Jan Jan 

where • denotes a placeholder for the argument over which the Green function is integrated. In 
order to obtain the representation formula for the vorticity, we remark that for any A E M^, 

V A(E-A) = VGA A, 

with G as defined above, so we obtain 


u}{x) = / 'VG{x — •) A / — / 'VG{x — •) A {T{u,p) —2u<S) ei) 
Jn Jan 

— / Va;(n • VG(a; — •)) A w + Va, (u • VG(a; — •)) A n 
The asymptotic expansions of the fundamental solutions are given by 


n 


B{x) = 


-1 / 1 


^/W2^ VrV2 


+ 0 


\ g—r(l—cos6) ( ^ 


., 3/2 


sm I 


+ 


1 (cos 0 sin ( 


VG(a;) = — (1 — cos 6, — sin0) 

V Stt 


U/2 


+ 0 


.. 3/2 


sin 6 1 — cos 6 I dvrr \ sin 9 — cos i 

—r(l—COS 0 ) 


and for |a| = 1, 


|D"E| < 


|D"VG| 


< 


( 1^1 ^ ^ -r(l-cos6») I 

i ^1/2 ^ j,3/2 J 

(\ 0 ? ^ ^ 


1/2 j.3/2 


r(l—COS 0 ) 


If / has compact support, the solution of the Oseen equation (2) behaves like the fundamental 
solution. 


“W = + 0(^375 + ;^) + O(^) ■ 

'V? , 1 ' 


x{x) = VG{x)-F^ + 0\ 


n/2 ' ^ 3/2 


r(l—COS 6 ) 


(3) 

(4) 
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where F is the net force 


F= f f+ f (T(-u,p) - 2u (g) ei) • n . 

Jn JdB 

Explicitly, the asymptotic expansion of the velocity is given by 

u{x) = U^{x) +Uh{x) g-r(l-cos0) + O 

where is the wake part and a harmonic function, 

^ -Eiei 1 

, . if cos 6 sin 0 \ ^ Fi Sr F 2 eg 

Uhlx) = -— . „ n \ F = --- . 

dvrr \ sm a — cos a I dvr r dvr r 


(5) 

( 6 ) 
(7) 


In contrast to the asymptotic expansion (3) or (5) of the velocity field, the asymptotic expansion 
(d) for the vorticity is only relevant inside the wake region, because outside the wake region, the 
remainder is larger than the asymptotic term. In order to obtain an asymptote that is relevant 
in all directions, we have to proceed differently, by using the following asymptotic expansion. 


VG{x - xq) = 

By applying this result, we obtain 

Uj{x) = VG{x) • F^{9) + o( — ^ ^-r(l-cos0) 


— r(l —cos 6 ) 


. 3/2 


( 8 ) 


where F{9) is now a function depending on the angle 9, 

F{9) = [ qMcos{9-9o)-cos9o)j^^ f -2ug)ei)-n. 

Jn JdB 

To our knowledge, this asymptotic formula is nowhere mentioned in the literature. With this 
expression, the asymptote is now detached in all directions from the remainder. 

3 Asymptote for the nonlinear problem 

For the nonlinear problem (1), neither the asymptotic behavior (5) nor (8) is correct, as shown 
below. The best results concerning the asymptotic behavior of u and a; for the nonlinear problem 
(1) are due to Babenko (1970). In particular he shows the following result: 

Theorem 1 (Babenko, 1970, Theorems 6.1 & 8.1). If u is a physically reasonable solution of 
(1), i.e. such that u — Uoo = for some e > 0 small enough, then the velocity satisfies 


u{x) — Uao = E{x)F + O 
and the vorticity satisfies 


1^1 (^^OS ^il) I —r(l—cos6) 


,1/2 +Vj'^''™”'+0vrl« 


ui{x) = VG(x) ■ + o(^) 

for some pL S (0,1), where E' E is the net force, 

F = I f+ I {T{u,p) - u^u) - n. 

Jn JdB 


d 










This result is optimal for the velocity in a sense that the remainder decays faster than the 
asymptotic terms. However, for the vorticity, this result is only relevant in the wake since the 
remainder is greater than the asymptotic term if 0 7 ^ 0, due to the fact that /r E (Ojl)- 
fact, we will show that the asymptotic behavior of the vorticity is not given by the Oseen tensor 
outside the wake. More precisely if we consider the vorticity which is given by 

Aoj — 2diOJ = u • Vw + V A / , 


then the Oseen approximation is given by it = 0. We will show that the linearization that leads 
to the correct asymptotic behavior of the nonlinear system is given by linearizing around the 
harmonic function of the Oseen tensor itself (7). We will show that this new linear system 
has an asymptotic behavior where the power of decay itself depends on F: 


Theorem 2. If there exists e E (0,1/2) and H, H E M sueh that the solution u of (1) satisfies 


\{u-Uoo- Uh) • er\ < -^e cosd )/2 


u 


^ 1/2 ' j,l+e 

r' /o r' 


\{u - - Uh) ■ ee\ < _e-Ki-cosO/2 + 


(9) 




where 


2A 2B 

Uh = 2V (Hlogr + BO) = —H-eg , 

r r 


for some v > 0 small enough, then the solution 0 / ( 1 ) with f a souree term of compact support 
satisfies for some C > 0, 


Iwl < (77-"4(l-cos6»)+Bsin6»^-l/2g-r(l-cos6») 


Moreover, 

^Jjl^x) = ^M^-cos9)+Bsme 
where p, is some 2TT-periodic funetion. 


rV 2 


+ 0 


nl/2+e 


r(l—cos 9) 


Remark 3. For Uoo, /, and u* small enough, the hypothesis (9) on u follows from theorem 1, 
with 


A = 


El 

Sir 


B = - 


Ei 

Stt 


Proof. Let u be a solution of (1) which by hypothesis can be written as u = u^o + Uh + u, with 
u satisfying 


\u-er\ < 


g-r(l-cos0)/2 ^ 


U+e ’ \U • eel _ • 

This expression is to be understood as To prove the result, we consider the vorticity equation 

Aw — 2diOJ — Uh • Vw = u • Vw + V A / . (10) 


< LQ-rL-<^ose)/2 


The change of variables, 

Uj{r,9) = ^^(l-cos 0 )+Bsineg-r(l-cose)^(^^^^ ^ 

transforms the original equation ( 10 ) into 

Aa — 2dra - = v • Va + (pa + R, (11) 

r 
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where v and (p are linearly related to Uh and u and satisfy 


\V • Br 


< 


n /2 ’ 


1 ^ • ee\ < 


z^logr 


and where the source term R is given by 

R(r,e) = r-^i^-^os9)-Bsin9^ril-cose) 


M< 




We will show that the solution o of (11) satisfies the bounds 


a 


< 

r\j 


,1/2 ’ 


Br • Va| < 


1 


Bg • Vo 



rs_/ 5 


r 


( 12 ) 


which imply the bound claimed on u. However in order to prove (12), a bootstrap argument in 
not sufficient. The idea is to consider the equation (11) as a linear equation in a with a and Vo 
given on dfl for fixed v, ip, and R and to construct a solution o by a fixed point argument. By 
uniqueness, the solution constructed by the fixed point argument is equal to the solution whose 
existence is assumed by hypothesis. 

The fundamental solution of the linear operator defining the left hand-side of (11) is 

W{x,xo) = ^e^~^°Ko{\x - xo\) , (13) 

Ztt 

where (ro, 0o) denotes the polar coordinates of xq. In view of the Green identity. 


J (^Ab + 2drb -\— J — 2c?j.a- ) ^ J — b^a + abB,.) • n , 

the solution can be written as 

a{x)= [ WR+ [ W{pa + vVa) + [ [aV- WVa + aWBr] ■ n , (14) 

Ju JQ JdQ 

where the integrations are performed over xq. For i = 1,2,3, we denote by ai{x) the ith term 
in the expression (14). 

The asymptotic expansion at large r, of the fundamental solution is given by 


+ O(^) , 


VxoW{x,xo) = 


(cos6>-cos6lo,sin(9-sin(9o) 1 




, 1/2 


+ 0 


,3/2 


and since R has compact support and dQ is a bounded, the first and the third term of (14) have 
the claimed asymptotic behavior. 


ai{x) 

asix) 


1 1 
^ri /2 


I 

Jm. 


,1/2 


+ 0 


^-ro[l-cos[9-9o)) ^ Q 




ri/2 



where pi and ps are 27r-periodic functions of the angle 9, 

^ f ^-M^-coseo)-Bsm9o^ro{cos{9-9o)-cos9o) ^ 

Jn ° 

[a (cos 0 — cos 00) sin0 — sin^o) — Vo-|-ae^] • n . 

vStt Jen 
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In the same way, 


V^W{x,xo) = 0 
and we obtain 

iVail < 


.3/2 


.3/2 ’ 


Va:Va;obh(®,®o) = O 


iVasI < 


.3/2 


.3/2 


SO oi and 03 satisfy (12). To deal with the second term 02 , we hrst make a hxed point on a in 
the space dehned by ( 12 ) so we have 


\ipa + vVal < 


j-S / 2-|-£i 


By using lemma 4 and lemma 5 we obtain that the second term 02 of (14) satishes (12). Since 
v is small, a hxed point argument shows the bound ( 12 ) claimed on a and therefore the bound 
on uj. 

In order to prove the asymptotic behavior of a 2 , we show that 


/ = 



,1/2 


,j-ro(l-cos(0-6lo)) 


3/2+e 


da^o 


is bounded by r By using the asymptotic expansion of the Bessel function Kq, we have 

/ < /i + I 2 where 


h = 
h = 


. 1/2 


^r-ro-ri 


. 1/2 


^-ro(l-cos( 0 - 0 o)) 


3/2+e 


dajo , 


L rf ‘ 


^r-ro-ri _ 


3/2+e 


+ , 


with ri = |a 3 — xq\. In view of lemma 5, I 2 is bounded by r For vq > r/2, the term Ii is 
bounded by . For rg < r/2, we have 


^ 1/2 j. 1/2 


1 ko ~ 2r cos (0 — 0o) I 

y.jJ^/ 2 ^ 1/2 T T\ 


ro 


< 

1/2 

r j- 


gr-ro-ri _ g-ro(l-cos(6»-0o)) ^ ^-mfl-cosfe-en)) rg sin {6 9q) 

ri + r — ro cos {0 — 9q) 

< „-ro(l-cos(0-0o))/o_l++(^ 

* 

r 

Therefore, always for ro < r/2, we have like in the proof of lemma 5, 


h< 


f 


^ 1/2 j.1/2 


ro(l—cos(0—0o)) _ 


1 + 




^r-ro-ri _ g-ro(l-cos(0-0o)) 


1 + r^/'+" 


X 2 

1 f ^ -ro(l-cos(e-en)) f + ^0 1 ^ ~ ^o| < 

r Lo 1/2 1 , 1/2+e -^0 ^ 

^ aR 2 1 + 


1 


.1/2+e 


□ 
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Lemma 4. The Green function (13) satisfies 


|VL| 


< 1 r-ri-rp 

^ 1/2 


|e^-VPL| 


'1 '1 
where ri = |£c — a;o|. 

Proof. First, the Bessel functions satisfy 


< _}_p(r-ri-ro)/2 
3/2^ ' 


1/2 

r-^ 


Ko{r) 


1 


Kfir) 


< 


1 


1 


Kfir) - Koir) < 


^1/2 ’ ^ ~ ^1/2 ’ 

SO the first bound is proven. Since r — tq cos(0 — 0o) = ^i cos(0i — 0), we have, 


avF = 


1 


27rri 


,,r-ro 


[(ro cos (0 - ffo) - r) Kfiri) + riKo(ri)] , 


^r—ro 


27r 


ri — r + ro cos{9 — 6o) 


ri 


Koiri) - cos( 6 » - 01 ) (i^i(ri) - Ko{ri)) 


so 


\drW\ < + -e^-'-a (|ro + ri - r| + tq |1 - cos(0 - 0o)|) Ko{ri) 

Ti' ^1 


1 


1 


< [1 + 1^^ + _ ^1 + |l _ cos(0 - 0 o)|] < — 


3/2 


3/2 


For the last bound, we have 


-deW = TT^e^ ^°rosin(0 - 0o)-^i(n) 
r ivrri 


so 




< _ , 

~ 3/2 


^r-ri-ro 


ro |sin (0 - 0 o)| < 


1/2 

< .(o_p(r’-ri-r-o)/2 


Lemma 5. For a, cr G (0,2) such that a + a > 3/2, we /lare 


ri r-Q 


1 


(j—3/2 


|logr| |logr| 


where r = |ai|, ro = |®o|, ri = |® — aiol CLnd 5 ^, 0 - denotes the Kronecker delta. 
Proof. We have to estimate 

r*oo r+TT 


I = 


roo r-\-7r ^ J 

JO J-TT ^0 


^r—ri—ro 


rod0odro . 


□ 


First of all, since 

r^ = r^ + rQ — 2 rro cos (0 — 0 o) > + ?’o cos ^(0 — 0 o) — 2 rro cos (0 — 0 o) = {ir — r^ cos (0 — 0 o))^ , 

we have 

?’i + J’o - > "ro (1 - cos (0 - 9q)) , 























and therefore 


/ +7r 




/. 


g-ro(l-cos(0-0o))(i0^ < 1 

— TT 


1 + 4'^' 


If ro < r/2, then ri > rj2 and therefore 

^ f-rl2 r+-K 2 


1 fT' ^ /■+7r 1 

I<- / -e-^°(i-“^(^-^“))rod0odro 

r“ io i-TT ^0 


' 0 j —TT ' 0 

1 1 1 


-f 

Jo 


0 1 + Tq' 


dr-o < 


|log r\' 


t,3/2 


1 


oOL^G — Zj^ 


In the case where ri < r j2^ we have by symmetry the previous bound with a and a exchanged, 


/ < 


|log r| 


1,3/2 


1 


rpG rpCX.~\~G — 3/2 

Therefore, it remains the case where tq > r/2 and ri > r/2, for which we have 


I < 


I 


1 


r/2 {r + |r-r„|)%J-V2 


^ dro < 




□ 
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